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Abstract 

We evaluate the ratio of the scale parameter A in domain wall QCD to the one in the continuum 
theory at one loop level incorporating the effect of massless quarks. We show that the Pauli- 
Villars regulator is required to subtract the unphysical massive fermion modes which emerge in the 
fermion loop contributions to the gluon self energy. Detailed results are presented as a function of 
the domain wall height M. 
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I. INTRODUCTION 



The domain- wall fermion formalism provides us an opportunity to realize the chiral sym- 
metry at finite lattice spacing with any number of quark flavors^, Q]. Since it is well known 
that the chiral symmetry plays an essential role in low energy QCD, this point is a superior 
feature over the Wilson and the Kogut-Susskind quark actions to perform lattice QCD sim- 
ulations. While in the past two decade both conventional quark actions have been widely 
used to make unquenched QCD simulations incorporating the dynamical quark effects, re- 
cent progress in computational technologies allows us to embark on a full QCD simulation 
with the domain- wall fermion jjj. 

One primary step in the full QCD simulation is to determine the renormalized coupling 
constant, which is a fundamental parameter in QCD. At the one-loop level the renormalized 
coupling in the continuum MS scheme gjj^(fi) at the scale /i and that on the lattice go are 
related by 

11/ 22 \ / 4 \ 

S^y = 7 + v 9 + w log(H ) + Nf \ f " 4^ log(H J (1) 

where d g and df are renormalization constants: the former depends on the gauge action 
and the latter on the quark action. These constant terms determine the ratio of the lattice 
A parameter and the continuum one at the one-loop level. From the expressions of the A 
parameters, 



A/ 



>S"*«exp (-^) , (2) 



with 



^ = ^9h(p)r^^{-^ &j ) (3) 



* = is? (" - l N ') • (4) 



we obtain the ratio 



^■MS 



(a^)exp ( -— { iA - -2 } ) + O(g ), (6) 

MS ' 

exp(-^{d g + N f d f }\+O(g 2 ). (7) 
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Thus the finite renormalization between the coupling constants defined by different regu- 
larization schemes is to relate their A parameters. In Table HI we list the renormalization 
constant d gi which has been evaluated for various gauge actions DDDDQ 

In this paper we make a perturbative determination of df for the domain- wall fermion, 
which is a necessary ingredient to extract the renormalized coupling constant from the 
numerical simulation with domain- wall QCD. While t he p erturbative techniques are well 



developed in our previous publications 



BE 
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12. 



13. 



141 ]. a characteristic feature in this 



calculation is to require the Pauli-Villars field to subtract the unphysical massive fermion 
modes, which emerge from the fermion loop diagrams. 

This paper is organized as follows. In Sec. II we introduce the domain-wall fermion 
action with the Pauli-Villars regulator. We give their Feynman rules relevant for the present 
calculation. In Sec. Ill we evaluate the quark loop contributions to the gluon self energy in 
the continuum and on the lattice. The ratio of continuum and lattice A parameters is given 
in Sec. IV. Our conclusions are summarized in Sec. V. 

The physical quantities are expressed in lattice units and the lattice spacing a is sup- 
pressed unless necessary. We take SU(A^ C ) gauge group with the gauge coupling constant 
9- 

II. ACTION AND FEYNMAN RULES 



We explain the domain-wall fermion action accompanied with its Pauli-Villars regulator 
in the present calculation. For the fermion action we take Shamir's form l|: 



DW 



EE 

n s=lL" £t 

+^ (^(n) s (l + j 5 )i/;(n) s+1 + tf(n) s (l - 75 )^(n) s _i) + (M — 1 + 4r)^(n) s ^(n), 



+m^2 {ip{n) Ns P R il){n)x + ^(n)iP L ip(n) Ns 



where is the link variable of the SU(A' r c ) gauge group. We label a four-dimensional space- 
time coordinate by n and an extra fifth- dimensional one by s, which runs from 1 to N s . Since 
no gauge interaction is imposed along the fifth dimension, we are also allowed to consider 
that the index s labels the N s "flavor" space. The parameter m denotes the physical quark 
mass and at m = one chiral zero mode is realized under the condition < M < 2 for the 
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Dirac "mass" (domain-wall height) M. The Wilson parameter is set to r 



-1. P 



R/L 



arc 



projection operators defined by Pr/l — (1 ± 75)/2- 

The "physical" quark fields on the four-dimensional space-time are constructed in terms 
of the fermion fields at the boundaries, 



q(n) = P R ip(n)i + P L ij(n) Ns , 
q(n) = i)(n) Na P R + Tp(n)iP L , 



(9) 
(10) 



with which we can extract the zero mode of domain-wall fermion at m = 0. 

We also introduce the Pauli-Villars field to subtract the massive modes of the domain- 
wall fermion, which emerge from the quark loop diagram in proportion to N s in the large 
N s limit. The action is given by 



S 



PV 



<Sbw(^ -> 4>, i> -> 4>, m = -1), 



(ii) 



where <ft is the iV s -flavor Wilson-Dirac boson, satisfying the anti-periodic boundary condition 
in the fifth direction because of the special value m = —1 15] . Here we have to remark 



one point: no Pauli-Villars field has been introduced in the previous calculations for the 

12L I13L Il4j , since it requires no 



renormalization factors of the composite operators jil llfl 
quark loop diagram. 

For the gauge action we employ the following form: 
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Sgluon = — S C TlU P l + Cl H ^Urtg + C 2 ^ Tr ^chr + C 3 ^ Trf/ p's 

^ I plaquette rectangle chair parallelogram 



(12) 



where the first term represents the standard plaquette action, and the remaining terms are 
six-link loops formed by a 1 x 2 rectangle, a bent 1x2 rectangle (chair) and a 3-dimensional 
parallelogram. The coefficients Co, • • • , C3 satisfy the normalization condition 



c + 8ci + 16c 2 + 8c 3 = I. 



(13) 



There are several choices for these parameters: c\ = c 2 = C3 = (Plaquette), C\ = —1/12, 
c 2 = c 3 = O(Symanzik)0, [l(| ci = —0.331, c 2 = c 3 = O(Iwasaki), c\ = — 0.27, c 2 + c 3 = 
-0.04(Iwasaki') JlJ, ci = -0.252, c 2 + c 3 = -0.17(Wilson)18] and c x = -1.40686, c 2 = 
c 3 = 0(DBW2)[19J. The last four cases are called the RG-improved gauge action, whose 
parameters are chosen to be the values suggested by some approximate renormalization 



group analyses. Some of these actions are now getting widely used to simulate the domain- 



wall QCdQ l2l| . since they are expected to realize smooth gauge field fluctuations which 
approximate the continuum gauge field better than the naive plaquette action. 

To prepare the Feynman rule we develop weak-coupling perturbation theory by writing 
the link variable in terms of the gauge potential 

1 



U^n) = exp I igaT A A A 



n 



-J 1 



(14) 



where T A (A — 1, ... , N% — 1) is a generator of color SU(iV c ). Here we omit the expression 
of the gluon propagator and vertices, since they are not relevant for the present calculation. 
Quark-gluon vertices are identical to those in the N s flavor Wilson fermion. We need one- 
and two-gluon vertices for the present calculation: 



VxMp)st = -gT A {i 1 , 1 cos{qj2 + pj2)+rsm{qj2+pj2)}5 sU 



-gT A v 1 ^(q,p)S. 



st ■ 



(15) 
(16) 



V 2 A *(q,p)st = ^g^-{T A ,T B }{i 1 ^(qj2 + pj2)-r C os(qj2+pj2)}5^ sU (17) 



2 U 2 

11 ,2 1 frpA rpB 



-g'-{T A } T B }v 2 » u (q,P)Ssu 



(18) 



where p and q represent incoming and outgoing quark momentum, respectively. 

The fermion propagator is obtained by inverting the domain-wall Dirac operator. In 
the present one-loop calculation we will take iV s — >• oo to neglect the contribution of the 
0(e~ aNs ) terms with the positive constant a given in eq. (|2fi|) . To this end it is sufficient to 
consider the following component of the fermion propagator: 



= (^-i^ fM smp fi + W~) su G R (u,t)P + 
-i 1/1 smp fl + W + ) G L (u,t)P- 



(19) 
(20) 



where the sum over the same index is taken implicitly. and Gr/l are given by 

f-W{p) 1 \ 
-W(p) "-. 

1 

V m —W(p) ) 



W + (p) 



st 



(21) 
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G L {s,t) 



cosh(a) 

A 
F 



W~{p) 



f-w( P ) 

1 -W{p) 



St 



V i 

W(p) = 1 - M - r 2(1 - cosp M ). 



m \ 



(22) 



(23) 



Gh(M) = |[-(i- m 2 )(i-iy e -j e 



a\ n(-2N s +s+t) 



-2W sinh(a)m ( e «(-".+»-*) + e^-^-^) 



A 
F 



1 - m 2 ) (1 - W^e a ) e ~ a{s+t) 
+ Ae -*\s-t\^ ( 24 ) 

2> 



1 - m 2 ) (1 - lfe a ) e a( - 27Vs+s+ '- 2) - (1 - m 2 ) (l - We~ a ) e a ^ s - t+2) 



-2Wsinh(a)m ( e °(-N.+«-*) + e ^-^-s+t) 

l + W^ + E^sin 2 ^ 
2\W\ 



+ Ae - a \ s -t\ 



2W/sinh(a)' 

1 - e a W - m 2 (l - We" 



(25) 
(26) 

(27) 
(28) 



In case that W becomes negative the fermion propagator is given with the replacement 



e ±a . In the present calculation we consider only m = case. The Feynman rule 



for the Pauli-Villars regulator is obtained by setting the quark mass to m = — 1. 



III. FERMION LOOP CONTRIBUTION TO THE GLUON SELF-ENERGY 

We first present the calculation of the vacuum polarization diagram in Fig. [I] (a) employing 
a continuum regularization scheme with the massless quarks. Here we choose the naive 
dimensional regularization (NDR) as the ultraviolet regularization scheme. Both the Dirac 
matrices and the fermion loop momenta of the Feynman integrals are defined in D = 4 — 2e 
(e > 0) in this scheme. To regularize the infrared divergence we take the external gluon 
momentum p to be off-shell p 2 ^ 0. 

The gluon two point function in Fig. ^ (a) is expressed as 



After some algebra we obtain 



(29) 



K B (P) = - N fji^ S M^P 2 - PMftO + l M/^/V) + Acont} (30) 



with 



A 



cont 
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(31) 



where 1/e = 1/e — 7 + ln(47r) and /1 is the renormalization scale. In the MS scheme, the 
pole term 1/e should be eliminated. 

Let us turn to the calculation of the one loop corrections to the gluon self energy on 
the lattice. We illustrate the relevant diagrams in Fig. ^ (a), (b) for the contribution of 
the domain wall fermion and (c), (d) for the Pauli-Villars regulators. Each loop diagram is 
expressed by 



n£f b (p) = -JV> 



n£ AS (p) = -n 



n (2tt) 4 
- d 4 k 



n£f B (p) = +n 



1 1 (2vr) 4 
- d A k 



U (2tt) 4 
ll {d > (v) = +N f 



Tr 
Tr 
Tr 
Tr 



VA(k,p + k) uv S F (p + k) vs V£(p + k,k) st S F (k) tu \ , (32) 



V 2 A *(k,k) st S F (k) ts 



(33) 



v i»( k >P + k) uv S F {p + k) vs vf u (p + k, k) st S F (k) tv \ , (34) 



V 2 ^{k,k) st S F (k) ts 



(35) 



A A 

where V lM , an d 'S'b denote the one-, two-gluon vertices and the propagator for the 

Pauli-Villars regulators. The sum of these contributions is expressed by IT^f (p): 



n£?(p) = ^ ab {p) + ngr (p) + n^(p) + n^(p) 



\ AB , 



\ AB , 



d 4 £; 



tt (2tt) 



(36) 
(37) 



Since the function IP^f (p) has an infrared divergence for p 2 — > 0, we extract it by employing 



an analytically integrable expression I^(p) which has the same infrared behavior as J~ (p), 



t(27t; 



4 F 



- d 4 fc 



{li„ B (k,p)-I* B (k,p)} 



(38) 



where 



I* B (k,p) = 6(A 2 -k 2 )Tr 



1 



1 



(39) 



with a cut-off A (< tt). Now we find that the second term in the right-hand side of eq. 
is regular in terms of p, we can make a Taylor expansion, 



iC(p) 



— 7r (2tt) 



+ 
+ 
+ 



- d 4 k 



{^(M)-/;J?(M)} 



(2tt) 4 



(40) 



(2vr) 4 2^ A %a PA 
where p, X = 1, 2, 3, 4. 

Let us consider each term in eq. (|40j) . We first point out that the linear term in p is 
not allowed from the Lorentz symmetry. The integration of I^(k,p) — l^(k,0) is easily 
performed analytically, 

d 4 k 



,(2.y {WM-WM)} 



( W - PmP.) { 3 log(A 2 /p 2 ) + g } + WV3 



(41) 



As for n^f(0) = Iwf(k, 0), we find that this function vanishes. This can be shown as 
follows. For the fermion contribution, 

nf» (0) = +iV ; |5 A B _/ ^ T^Tr fc^^M*. fc)S' F (A:) rt ] , (42) 



2 



(27T) 

- d 4 k 



+N 



fTT°AB 



(2tt) 4 

- d 4 k 



Tr 
Tr 
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Uiit(A;, k)—S F (k) tt 
' d 

v ltl (k,k)S F (k) tt 



dk v 



-nS^(o), 



where we use the identities 



g 

g 

g 
dk y 



S F 1 {k) st = v lv (k,k)S : 



st • 



[S F l (k) st S F (k) tu ] = 0, 
v 1/x (k,k) = v 2 ^{k,k). 



In a similar way we can show 



4f B (o) = -ng>) 



(43) 

(44) 
(45) 

(46) 
(47) 
(48) 

(49) 



for the contribution of the Pauli-Villars regulator. After all, the sum of all the contributions 
vanishes: 

iOo) = n£f >) + ng A >) + m A >) + n£f >) = o. (so) 
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We now turn to the last term in eq. (pfUj) . The calculation of this term requires tedious 
algebra and yields quite lengthy results that make it practically difficult to give their full 
expression here. However, we think that it would be instructive to explain why the Pauli- 
Villars regulator is required by showing a part of the calculation. Let us take the term 
proportional to Tr [7^70,7^7^] as an example, 

ffl ( k 'P) = -^/#o^ Tr [7M7a7^7/3]cos(/^ + p ll /2)sm(k a + p a )cos(k u + p u /2)sm(k f3 ) 



x {G' R (s, t)G R (t, s) - G' R (s, t)G R {t, s)} 



(51) 



where G' R and G R are functions of p + k, while Gr and Gr are functions of k. Gr is obtained 
from eq. (|24J) with m = — 1. The fermion contribution is evaluated as 



G' R (s, t)G R (t, s) = £ [^e- a ' |s -" + A' L e a ' is+t - 2N ^ + A' R e a '^ s - t+2) 



s,t 
X 



[ Ae ~a| S -*| + ALe a(s+t-2N s ) + 

( —O p -{ a '+ a ) 1 _|_ p -(a'+a) ) 

A'Al _ . _ , 2 +iv/ 1 



(52) 



[1 _ e -(a'+a)] 

+A'(A L + A R ) 



' 1 _ e -(a'+Q!) 



1 + e (Q ' +a) e 2a 



Q 2a 



+(A' L + A' R )A 



e {a'+a) I 
1 + e (a'+a) e 2a' 
1 _ e (a'+a) 1 _ e 2a' 



+(A' L A L + A' R A R ) 



3 2(a'+a) 



[1 _ e (a'+a)]' 



with 



1 - We~ a 
1 - We a 

-e- 2a A, 



A. 



(53) 

(54) 
(55) 



where the superscript / denotes the function of p + k. On the other hand, the Pauli-Villars 
part is given by 



Gf R {s,t)G R {t,s) = 



e -a'\s-t\ _ e a'(s-t-N s ) _ e a'(-s+t-N a ) 



XA j^ e - Q l s -* _ e <x(s-t-N 3 ) _ e a(-s+t-N s ) 



A'A 



_2 e ~( a '+ a ) 1 -l e -( a '+ a ) 2e~( a ' +a ^ 
+ N s - -— T + 



(56) 
(57) 



^ [1 _ e -(ai+a)y ■ "l_ e -(a'+a) ' ^ _ g-( a '+a)] 2 j ' 

where we find the exactly same 0(N S ) term as in the domain- wall fermion part. This is 
sustained even after taking the derivative in terms of p according to eq. ()40|) . Through this 



kind of analytical calculation we have checked that all the 0(N S ) terms of the fermion 
contribution in the last term of eq. PUJ) are canceled out by the Pauli-Villars one, which is 
proportional to N s after collecting all the terms. 
After all, eq. (|4(J|) is summarized as 

r 2 \ cr^l 



+A5^p 2 + Bp^p, 



7T 



§ log I > . 

3 \p z a z 



(58) 



The momentum integration is carried out numerically employing the Monte Carlo routine 



BASES 22J. We choose A = ir for the cut-off. The results are checked by a mode sum for a 
periodic box of a size L 4 with L = 128 after transforming the momentum variable through 
k'^ = — sinfc^. The estimated values for A + 1/3 and B are presented in Table ITT1 as a 
function of M. We observe that the expected relation A + 1/3 = —B is well satisfied within 
error bars. Finally we obtain 



Ui B (p) = -N 



where 



(47T) 



A 



VvVv 



--log(pV) 



A 



latt 



la,;t = A + - + - l0g(7T 2 ) 



(59) 



(60) 



Comparing the above expression to the continuum counterpart in eq.(|30j). we obtain the 
finite renormalization constant df in eq.(|T|). 

1 



d t 



167T^ 



(Ai att - A 

cont ; 



(61) 



where the numerical values are given in Table |H] as a function of M. We find that the 
sign is opposite to and the magnitude is compatible with the Wilson, the clover and the 



Kogut-Susskind cases : df 
respectively. 



0.0066949 



d f = 0.0314917 



, Uf 



and d f = 0.0026248 



IV. THE RATIO OF CONTINUUM AND LATTICE A PARAMETERS 

We first show how to obtain the mean-field improved MS coupling g^(fi) at the scale \i 
from the lattice bare coupling g\. A detailed description for the quenched case is already 
given in Sec. VI of Ref. [lj]. For the full QCD case the mean field improvement modifies 
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the relation of and g$ in eq. (JTJ) : 



IP 22 / 4 

-2 + d g + c p + -— log(//o) + JV> ( d/ - —-3 log(/ia) ) . (62) 



(^(A*) </§ ^ ^ p 16^ ^ ' > 48vr 
One may use an alternative formula 27 1 

1 c P + 8c x Pl + 16c 2 P2 + 8c 3 P3 



22 

+ (C • C p + 8Ci • Cfli + 16C 2 ■ C R2 + 8c 3 • C R3 ) + — -r log(//d) 

l07T 

+^/(rf/- i ^log(/ia)), (63) 



where 



P = ^ Trf/ P Za ^ette = 1 - c P 5-o + O(^o)> ( 64 ) 
PI = -TrC/ rectang i e = 1 - c m #o + O(#o), (65) 
P2 = -TrUchair = 1 - c R2 gl + O(# 4 ), (66) 

P3 = T^vU parallelogram = 1 ~ CR 3 #q + 0(# 4 ), (67) 

and the measured values may be employed for P, PI, P2 and P3. Table [Q summarizes the 



values o 
of Ref. 



Sj>i c Rii C R2 an d C R3 f° r various gauge actions, which are taken from Table XVI 



14 1 . One may also use the mean-field estimate for M, which is given by M 



M — 4(1 — P 1 / 4 ) 14|. No other modifications are necessary in this case: Just use df at M. 

As discussed in Sec. HJ we can determine the ratio of and A-^ from the finite renor- 
malization of the coupling constant. From eq.((7j), the ratio is written as 

where the first factor is given in Table U for the various gauge actions, and the second one 
is given in Table ITT1 as a function of M, in the case that Nf = 0, 1, 2, 3, 4. 



V. CONCLUSION 



In this paper we evaluate the fermion contribution to the finite renormalization of the 
coupling constant in the domain- wall QCD, which determines the ratio of the lattice A 
parameter to the continuum one. The numerical values are given as a function of M. We 
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show that the massive fermion modes emerging from the fermion loop contribution to the 
gluon self energy can be subtracted by the Pauli-Villars regulator. We explain how to 
implement the mean field improvement for the renormalized coupling <?fjg(/-0 and the A 
parameter ratio. 
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TABLE I: Values of the perturbative quantities d g , cp, and cm for the various gauge actions. Gluon contribution to the A parameter ratio 
in eq,(|68[) is given for the case of Nf = 0, 1, 2, 3, 4. 



action 




C3 


dg C p CR1 CR2 CR3 


N f = 


JV> = 1 


exp(-V(2/?o)) 
Nf = 2 Nf = 3 


N f = 4 


Plaquette 








-0.4682 1/3 0.5747 0.5228 0.5687 


28.81(2) 


35.78(3) 


45.80(4) 


60.80(5) 


84.45(8) 


Iwasaki 


-0.331 





0.1053 0.1401 0.2689 0.2223 0.2465 


0.4696(3) 


0.4473(3) 


0.4231(3) 


0.3970(3) 


0.3687(3) 


DBW2 


-1.40686 


0.5317 0.0511 0.1040 0.0818 0.0921 


0.02200(2) 


0.01720(1 


) 0.01300(1) 


0.009423(8) 


0.006488(6) 


Symanzik 


-1/12 





-0.2361 0.2440 0.4417 0.3846 0.4215 


5.445(4) 


6.074(5) 


6.879(6) 


7.935(7) 


9.365(9) 


Iwasaki' 


-0.27 


-0.04 


0.1471 0.2797 0.2342 0.2611 












Wilson 


-0.252 


-0.17 


0.1196 0.1286 0.2439 0.2070 0.2352 


0.4238(3) 


0.4010(3) 


0.3765(3) 


0.3502(3) 


0.3220(3) 



TABLE II: Numerical values for A + 1/3, B and df as a function of M Fermion contribution to 
the A parameter ratio in ea. ()68|) is given for the case of Nf = 1, 2, 3, 4. 



A + 1/3 



B 



d 



L. 



eM-Nfd f /(2P )) 



N 



L 



1 



N 



L 



N 



N 



L 



-5.4253(27) 

-4.6075(22) 

-4.1408(20) 

-3.8102(18) 

-3.5497(18) 

-3.3341(15) 

-3.1475(15) 

-2.9819(14) 

-2.8310(13) 

-2.6925(12) 

-2.5635(12) 

-2.4418(12) 

-2.3272(11) 

-2.2176(10) 

-2.11237(99) 

-2.01129(98) 

-1.91395(88) 

-1.82014(87) 

-1.73060(81) 

-1.64359(81) 

-1.56012(77) 

-1.48027(71) 

-1.40379(68) 

-1.33232(65) 

-1.26633(62) 

-1.20519(59) 

-1.15142(57) 

-1.10734(55) 

-1.07434(53) 

-1.05539(52) 

-1.05344(53) 

-1.07444(54) 

-1.12617(56) 

-1.21873(61) 

-1.36869(68) 

-1.60462(80) 

-1.97527(99) 

-2.5818(13) 

-3.7162(19) 



5.4259(27) 

4.6092(23) 

4.1414(19) 

3.8106(18) 

3.5494(18) 

3.3346(16) 

3.1483(14) 

2.9825(13) 

2.8316(12) 

2.6938(12) 

2.5650(12) 

2.4433(11) 

2.3282(11) 

2.2182(11) 

2.11208(98) 

2.01132(97) 

1.91450(89) 

1.82057(91) 

1.73063(86) 

1.64365(82) 

1.55976(76) 

1.47976(71) 

1.40313(69) 

1.33289(66) 

1.26646(62) 

1.20600(60) 

1.15256(57) 

1.10903(55) 

1.07514(54) 

1.05533(52) 

1.05437(53) 

1.07497(54) 

1.12665(56) 

1.21856(61) 

1.36952(68) 

1.60601(80) 

1.97602(99) 

2.5816(13) 

3.7153(19) 



-0.028209(17) 
-0.023030(14) 
-0.020075(13) 
-0.017981(11) 
-0.016332(11) 
-0.0149658(97 
-0.0137846(96 
-0.0127360(87 
-0.0117805(83 
-0.0109029(79 
-0.0100864(76 
-0.0093153(73 
-0.0085901(67 
-0.0078955(64 
-0.0072294(63 
-0.0065893(62 
-0.0059729(55 
-0.0053788(55 
-0.0048118(51 
-0.0042608(51 
-0.0037322(49 
-0.0032266(45 
-0.0027422(43 
-0.0022896(41 
-0.0018718(39 
-0.0014846(37 
-0.0011441(36 
-0.0008650(35 
-0.0006560(34 
-0.0005360(33 
-0.0005236(33 
-0.0006566(34 
-0.0009842(36 
-0.0015704(39 
-0.0025200(43 
-0.0040140(51 
-0.0063612(63 
-0.0102019(82 
-0.017386(12) 



1.24054(16) 
1.19240(13) 
1.16578(11) 
1.14728(10) 
1.132911(97 
1.121148(83 
1.111075(81 
1.102208(73 
1.094190(69 
1.086878(65 
1.080117(62 
1.073772(60 
1.067839(55 
1.062186(52 
1.056794(51 
1.051638(50 
1.046696(44 
1.041956(44 
1.037451(40 
1.033092(41 
1.028928(39 
1.024961(35 
1.021174(33 
1.017649(32 
1.014405(30 
1.011408(29 
1.008780(28 
1.006631(27 
1.005025(26 
1.004104(25 
1.004009(26 
1.005030(26 
1.007548(27 
1.012071(30 
1.019442(34 
1.031146(40 
1.049807(50 
1.081071(68 
1.14207(10) 



1.58537(44) 
1.45676(33) 
1.38810(29) 
1.34143(25) 
1.30577(24) 
1.27696(20) 
1.25255(20) 
1.23128(17) 
1.21221(16) 
1.19496(15) 
1.17912(15) 
1.16436(14) 
1.15064(13) 
1.13767(12) 
1.12535(12) 
1.11365(11) 
1.10249(10) 
1.091844(98 
1.081777(90 
1.072083(90 
1.062866(85 
1.054123(78 
1.045815(73 
1.038112(70 
1.031050(66 
1.024549(63 
1.018865(60 
1.014230(58 
1.010774(56 
1.008794(55 
1.008590(55 
1.010784(56 
1.016207(59 
1.025985(65 
1.042025(74 
1.067771(89 
1.10951(11) 
1.18135(16) 
1.32845(26) 



2.10102(95 
1.83331(68 
1.69612(57 
1.60520(48 
1.53700(45 
1.48273(38 
1.43734(36 
1.39822(32 
1.36349(30 
1.33236(28 
1.30403(26 
1.27784(25 
1.25368(22 
1.23097(21 
1.20958(20 
1.18937(19 
1.17023(17 
1.15207(17 
1.13501(15 
1.11867(15 
1.10321(14 
1.08863(13 
1.07484(12 
1.06211(12 
1.05050(11 
1.03985(10 
1.030569(98) 
1.023026(94) 
1.017416(90) 
1.014206(89) 
1.013876(89) 
1.017431(91) 
1.026241(96) 
1.04220(11) 
1.06857(12) 
1.11143(15) 
1.18225(19) 
1.30800(28) 
1.58024(49) 



2.9127(19) 
2.3937(13) 
2.1400(10) 
1.97679 
1.85699 
1.76331 
1.68612 
1.62042 
1.56279 
1.51167 
1.46560 
1.42339 
1.38480 
1.34883 
1.31520 
1.28368 
1.25404 
1.22612 
1.20005 
1.17525 
1.15194 
1.13008 
1.10952 
1.09065 
1.07352 
1.05788 
1.04431 
1.03332 
1.02517 
1.02052 
1.02004 
1.02520 
1.03800 
1.06132 
1.10022 
1.16431 
1.27263 
1.47203 
1.93268 
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FIG. 1: One loop corrections to the gluon self energy. Solid line is for the propagator of the 
domain- wall fermion and broken one for the Pauli-Villars regulator. 
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